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Temperature anisotropies in the Cosmic Microwave Background (CMB) are affected by the late 
Integrated Sachs- Wolfe (IISW) effect caused by any time-variation of the gravitational potential on 
linear scales. Dark energy is not the only source of IISW, since massive neutrinos induce a small 
decay of the potential on small scales during both matter and dark energy domination. In this work, 
we study the prospect of using the cross-correlation between CMB and galaxy density maps as a tool 
for constraining the neutrino mass. On the one hand massive neutrinos reduce the cross-correlation 
spectrum because free-streaming slows down structure formation; on the other hand, they enhance 
it through their change in the effective linear growth. We show that in the observable range of 
scales and redshifts, the first effect dominates, but the second one is not negligible. We carry out 
an error forecast analysis by fitting some mock data inspired by the Planck satellite. Dark Energy 
Survey (DES) and Large Synoptic Survey Telescope (LSST). The inclusion of the cross-correlation 
data from Planck and LSST increases the sensitivity to the neutrino mass by 38% (and to the 
dark energy equation of state w by 83%) with respect to Planck alone. The correlation between 
Planck and DES brings a far less significant improvement. This method is not potentially as good 
for detecting rui, as the measurement of galaxy, cluster or cosmic shear power spectra, but since it 
is independent and affected by different systematics, it remains potentially interesting if the total 
neutrino mass is of the order of 0.2 eV; if instead it is close to the lower bound from atmospheric 
oscillations, ~ 0.05 eV, we do not expect the ISW-galaxy correlation to be ever sensitive to m^. 

PACS numbers: 98.80.Cq 



> 
in 

in 
in 

d 
o 



I. INTRODUCTION 



As photons pass through a changing gravitational po- 
tential well, they experience a redshift or a blueshift, de- 
pending on whether the well grows or decays respectively. 
Cosmic microwave background (CMB) photons can expe- 
rience such variations between the time of last scattering 
and their detection now. This effect was first described 
by Sachs and Wolfe in 1967 [H, and hence is dubbed 
the integrated Sachs- Wolfe effect (ISW). During a Cold 
Dark Matter (CDM) and/or baryon dominated era, the 
gravitational potential distribution remains frozen, and 
the ISW effect has no net effect on the blackbody tem- 
perature of CMB photons. This property is crucially re- 
lated to the fact that non-relativistic matter (like CDM 
and baryons) has a vanishing sound speed, and experi- 
ences gravitational clustering on all sub-Hubble scales af- 
ter photon decoupling, as described by the Poisson equa- 
tion. In such a situation, the universal expansion and the 
gravitational contraction compensate each other in such a 
way as to maintain a static gravitational potential. How- 
ever, when the expansion rate is affected by any type of 
matter with a non-vanishing sound speed, e.g. during 
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Dark Energy (DE) domination, the gravitational pertur- 
bations decay and the cosmic photon fluid experiences a 
blue shift, acquiring extra temperature perturbations re- 
lated to the intervening pattern of matter perturbations. 
It was first proposed by Crittenden and Turok in 1995 Q 
to cross correlate maps of temperature perturbations in 
the CMB with those of matter overdensities in large scale 
structures (LSS), in order to measure a possible accelera- 
tion of the universe's expansion. However, the CMB and 
LSS data available at that time were not good enough 
for such an ambitious goal, and the first strong indica- 
tion of a positive acceleration came in 1998 from the side 
of type-la supernovae [1,0]. Analyses of the first (2003) 
and second (2006) data releases of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) Q were the first 
to indicate the existence of Dark Energy independent of 
acceleration, by means of the location of the second peak 
in the CMB power spectrum. Simultaneously, a number 
of interesting papers presented the first detections of the 
ISW effect by cross-correlating WMAP aniso tropy m aps 
with various LSS data sets 0,1, i, 0, [O, E Ellli, Ei , 
now able to give an independent measure for the accel- 
eration of the expansion of the universe. 

The domination of Dark Energy is not the only source 
of gravitational potential evolution and of a net ISW 
effect. On small cosmological scales, as soon as mat- 
ter perturbations exceed the linear regime, gravitational 
perturbations start to grow and to redshift CMB pho- 
tons. This effect, called the Rees-Sciama effect, has 
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not been significantly detected until now [l^. CMB 
jhotons can also be scattered by gravitational lensing 
1 3 an d by the Sunyaev-Zeldovich (SZ) effect [H] (see 
br detections in CMB-LSS cross-correlation 
analysis). An other party expected to affect the evolu- 
tion of gravitational perturbations -at least by a small 
amount- is the background of massive neutrinos. Over 
thirty years ago massive neutrinos were proposed as a 
Hot Dark Matter (HDM) candidate, and later ruled out 
as the dominant dark component, since HDM tends to 
wash out small scale overdensities during structure for- 
mation [2^ . Observed neutrino oscillations however con- 
strain neutrinos to have a mass [21I, In addition, 
the presence of a Cosmic Neutrino Background (CNB) 
is strongly suggested on the one hand by the abundance 
of light elements produced during primordial nucleosyn- 
thesis [H, [2I [23, and on the other hand by CMB 
anisotropics [H,!!, [H, [H, M, El, [H, [H, [H, [H . There- 
fore, a small fraction of HDM is expected to coexist with 
the dominant CDM component. On small cosmological 
scales (for instance, cluster scale), the free-streaming of 
massive neutrinos should induce a slow decay of gravita- 
tional and matter perturbations (ssj , acting during both 
matter and Dark Energy domination. This effect depends 
on the total neutrino mass summed over all neutrino fam- 
ilies, TOjy = X^i"^*! unlike laboratory experiments based 
on tritium decay or neutrinoless double-beta decay, which 
probe different combinations: hence, a cosmological de- 
termination of the total neutrino mass would bring com- 
plementary information to the scheduled particle physics 
experiments [13, [11] . The free streaming of massive neu- 
trinos has not yet been detected [11], but there are good 
prospects to do so in the future, since the smallest total 
neutrino mass allowed by data on atmospheric neutrino 
oscillations {nii, > a/ Am^tm ^ 0.05 eV) implies at least 
a 5% suppression in the matter/gravitational small-scale 
power spectrum [stI [ssj . A positive detection -even in 
the case of minimal mass- could follow from the analy- 
sis of future galaxy /cluster redshift surveys [sol. [iol. W\ . 
weak lensing surveys [4^. [43j. Lyman-a forest analysis, 
cluster counts [4^], etc. The goal of measuring the neu- 
trino mass from cosmology is very ambitious since each 
of these methods suffers from its own source of system- 
atics (bias issues, modeling of non-linear clustering, ...). 
Therefore, a robust detection could only be achieved by 
comparing the results from various types of experiments. 

The goal of this work is to describe a possible cos- 
mological determination of the absolute neutrino mass 
scale through the ISW effect induced by neutrino 
free-streaming on CMB temperature maps, using as 
an observable the cross-correlation function of galaxy- 
temperature maps. This possibility was investigated pre- 
viously by Ichikawa and Takahashi [3] (and suggested 
again recently in [1^). As neutrinos slow down the 
growth of structure, we expect the blueshift caused by 
an accelerated expansion to be more pronounced if neu- 
trinos have a larger mass. On the other hand, the distri- 
bution of matter inducing the late ISW effect is smoother 



in case of free-streaming by massive neutrinos. These two 
antagonist effects should in principle induce some mass- 
dependent variations in the galaxy-temperature cross- 
correlation function. 

In section|lI]of this paper we give an outline of the the- 
ory of the ISW-cffect in the presence of a neutrino mass. 
In section IIIIl wc use some mock data with properties 
inspired from the Planck satellite, Dark Energy Survey 
(DES) and Large Synoptic Survey Telescope (LSST) in 
order to show the potential impact of this method in the 
future. 



II. THE GALAXY-ISW CORRELATION IN THE 
PRESENCE OF NEUTRINO MASS 

A. Definitions 

The observed galaxy overdensity 6g in a given direc- 
tion n is defined as 



Sain) = / dzb{z)(j)Giz)6,„{n,z), 



(1) 



where z denotes redshift, b{z) is the redshift dependent 
bias function relating the observed galaxy overdensity to 
the total matter overdensity, and 4>g{z) is the galaxy 
selection function which can be chosen such that only 
galaxies within a certain range of redshift are considered. 
The observed CMB temperature map 



T{n) - To 



(2) 



results from various contributions, classified as primary 
or secondary anisotropics. By definition, secondary 
anisotropics are induced after photon decoupling and can 
be correlated to some extent with the surrounding large 
scale structure. The ISW component is one of these 
terms, and can be obtained by integrating the scalar met- 
ric perturbations (or just the Newtonian gravitational 
potential on sub-Hubble scales) along each line-of-sight 
between the last scattering surface and the observer. If 
the gravitational potential is written as a function of di- 
rection h and redshift z, the ISW term reads 



(n) 



dz—{n,z) 
az 



(3) 



where Zdcc is the redshift at decouplig. Immediately after 
decoupling and before full matter domination, the grav- 
itational potential does vary with time: this is known 
as the early ISW (elSW) effect, in contrast with the 
late ISW (IISW) in which we are presently interested. 
The two maps can be computed sep- 

arately by cutting the above integral in two pieces at 
some intermediate redshift chosen during full matter 
domination, when the gravitational potential is static. 
Note that in presence of massive neutrinos, the poten- 
tial is never really static on small scales, so the quan- 
tity A^'^^ might not be uniquely defined. Anyway, 
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this question is not relevant in practice. The observable 
quantity is not the late ISW auto-correlation function 
{Aif^'^{n)Aif^'^{n')), but only its cross-correlation with 
a given survey {Aif^^{n)SG{fi')). Then, the redshift dis- 
tribution 4)g{z) selects the range in which the ISW effect 
is being probed, and the choice of z» becomes irrelevant 
provided that remains larger than the redshift of all 
objects in the survey: (/)g(z,) ~ 0. 

Assuming that the galaxy-temperature cross- 
correlation function arises solely from the late ISW 
effect (i.e., assuming that other secondary anisotropics 
potentially correlated with LSS can be separated or 
have a negligible amplitude, which is a good assumption 
on the scales considered hereafter), we can relate the 
galaxy-temperature correlation multipoles to the real- 
space correlation function [Al]^^^ {n)6G{n')) . In the 
Limber approximation (see Appendix), one gets 
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where r(z) is the conformal distance up to redshift 
z, i7o = lOO/i km/s/Mpc is the Hubble parameter 
today, and the matter power spectrum is defined as 
{5m{k,z)5m{k' ,z)) = P{k,z) S^{k- k'). Note that we 
used the Poisson equation in flat space in order to relate 
the gravitational potential $ to the matter overdensity 
6m, and assumed a(0) = 1 by convention. Finally, the 
multipoles Cf^^ define the angular correlation function 
in a Legendre polynomial basis (pi), 
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Eq.(|3]) is often written in a form which assumes that 
the matter power spectrum is a separable function of 
wavenumber and redshift. This applies to the case of 
a (flat) ACDM universe, for which one can write 



z) = L>(A; zfa{zfP{k, 0) 



(6) 



with dzD = during full matter domination and d^D > 
during A domination. Figure [T] (left) shows the evolu- 
tion of D as a function of z for fl^ = 1 — 17,„ = 0.69. 
In the case of time- varying Dark Energy, the situation is 
qualitatively similar, and D just depends on more free 
parameters than A. In the rest of this paper, we will just 
write this function as D{z) for concision. 



Cf^, it is possible to use some approximate solutions 
valid only on the largest and smallest wavelength (see [38j 
and [i^ for more details). First, for wavelengths larger 
than the maximum value of the neutrino free-streaming 
scale, reached at the time of the transition to the non- 
relativistic regime, the power spectrum P-^" is completely 
unaffected by neutrino masses, and identical to that in 
a massless neutrino model with the same cosmological 
parameters (in particular, the same Qm and h) noted as 

pO. 

Vfc < fc„„ P-^- (fc, z) = [D{z)a{z)]'^P^- (fc, 0) 

with P-^"(fc,0) = P°(fc,0) . (7) 

On the other hand, for wavelengths smaller than the the 
free-streaming scale today, both the linear growth factor 
and the amplitude today are affected by neutrino masses, 
approximately like: 

Vfc > fcfa, Pf-{k,z) ~ [D{z)a{z)]^-if'^pf"{k,0) 
with P-^'^ (fc, 0) ~ [1 - 8/^]P°(fc, 0) , (8) 

where fi, = fl^/fl„i stands for the neutrino density to- 
day relative to the total matter density (so f2„i includes 
baryons, hot and cold dark matter). Here D{z) is always 
the same function, computed either for f„ 7^ on large 
scales, or for = on any scale, with a common value 
of Q\ (or of Dark Energy parameters) . The first approx- 
imation in Eqs. ([8]) is very accurate, as shown in Fig. [T] 
(left) where we compare the precise linear growth factor 
obtained numerically with the above solution. The sec- 
ond approximation is poorer, but more accurate ones can 
be found e.g. in Refs. fssl. lisj. 

Assuming that the galaxy selection function is very 
peaked around a median redshift Zm, the multipole Cp^ 
probes mainly fluctuations around the scale fc ^ l/r{zm). 
If / is larger than fcfsr(zm), Cf^^ is affected by neutrino 
masses through the term between brackets in Eq. 
Using Eqs. (O, this term varies with like: 

d.^^^^^[il + C{z)U){l^&U)]d.^!^ , (9) 
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For a typical Dark Energy model, the density fraction 
i^DE becomes negligible for z > 2 [sH], and hence the 
ratio D'/D is tiny. So, at hight redshift, the net effect 
of the neutrino mass is to increase the integrand in Cf~^ 
like: 



B. Effect of neutrino masses 

In models with massive neutrinos, the spectrum is not 
a separable function anymore (in other words, the lin- 
ear growth factor is scale-dependent), and Eq.Q cannot 
be further simplified. However, in order to make an- 
alytical estimates of the impact of neutrino masses on 
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This just reflects the fact that at high redshift, the ISW 
effect would be null on all scales for fi, = 0, while for 
/i/ > it is still active on small scales. However, for 
z < 2, D'/D becomes larger, and for typical values of 
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FIG. 1: (Left) Redshift evolution of the small-scale linear 
growth factor, defined here as [P{k, z)/ P(k,0)]^^^ /a{z) for 
k ~ lO/tMpc"^, and obtained numerically with CAMB for 
— 0.69. The lower curve corresponds to — and is 
exactly equal to the quantity D{z) defined in Eq.©. The up- 
per, solid curve corresponds to fu = 0.1, and is well approx- 
imated by the dotted curve, which corresponds to the first 
of Eqs.®. (Right) The function C{z), defined in Eqs.dU, 
computed here for Qa ~ 0.69. Roughly speaking, the effect 
of neutrino masses on Cj^*^ changes of sign when this function 
crosses eighth. 



Hue ^ 0.7 there is always a redshift below which C{z) is 
smaller than eight. Then, the term between brackets in 
Eq. ^ is smaller than one, and the net effect of neutrino 
masses is to decrease dziP/a^]- In Fig. [T] (right), we plot 
the fimction C{z) in the case of a cosmological constant 
with ri\ = 0.69. We see that C ^ 8 for z ~ 2; so, 
around this redshift and for I > fcfsr(z„j), the net effect 
of neutrino masses on Cp"^ changes of sign. 

In summary, if Zm is small, the expected effect of 
neutrino masses on the cross-correlation multipoles Cp^ 
consists in a step- like suppression at large Ts, quali- 
tatively similar to that observed in the galaxy auto- 
correlation multipoles Cp*^. However the suppression 
factor is smaller, since the lack of power in the matter 
power spectrum caused by neutrino free-streaming is bal- 
anced by the excess of ISW effect due to the behavior of 
the linear growth factor in presence of massive neutrinos. 
When Zm increases, the boost related to the ISW effect 
is seen more clearly, and ultimately, when Zm is chosen 
before dark energy domination, the net effect of neutrino 
masses is to increase Cp^ at large I. 

In order to check and quantify these effects, we com- 
puted the cross-correlation multipoles Cf^^ (and also for 
comparison the auto-correlation multipoles Cp*^) for two 
different cosmological models, sharing the same param- 
eters fib = 0.053, n,n = 0.31, r^A = 0.69, h = 0.65, 
A = ln[10lOA:37^2]^^ool/Mpc = 3.16, = 0.95, but 
with two different values of the neutrino density frac- 
tion = fl^/^rn, equal to or 0.1 (corresponding to 
three neutrino species sharing the same mass = or 
nil, — 0.41 eV). We adopted a galaxy selection function 



FIG. 2: Ratio of the cross-correlation multipoles C; and 
auto-correlation multipoles Cp*^ obtained for two cosmologi- 
cal models with neutrino density fractions equal to = 0.1 
or 0, and the same value of other cosmological parameters 
(see the text for details). 



of the form 
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peaking near the median redshift Zm = IAzq. For il- 
lustrative purposes, we choose the four values Zm = 
0.1,1,2,3, although in practice it would be very chal- 
lenging to map Sain) for z > 2: presently, available data 
with a reasonable signal-to-noise ratio range only from 
z - 0.1 to z - 1.5. 

In Fig. [2] we plot the ratio of the multipoles C^^ in the 
two models, compared with the same ratio for Cp*^. The 
free-streaming of massive neutrinos is responsible for the 
step- like suppression of Cp*^, like in the power spectrum 
P{k). The value of z™ controls the angle under which the 
free-streaming scale is seen in the map Sain), and hence 
the scale at which the suppression occurs in multipole 
space. As expected from the previous discussion, the 
neutrino mass effect on Cf^'^ is similar to that on Cp*^ for 
small Zm < 1, although the suppression factor is slightly 
smaller, due to the excess of ISW effect in presence of 
massive neutrinos. For z^ > 1, the amplification effect 
due to this excess has a clear and distinct signature at 
I > 100, and for z,„ ~ 2 the ratio displayed in Fig. ([2]) 
has a dip around I ~ 150. Unfortunatly, we will see in 
Sec. Ill Ci that for Z > 100 this effect is masked by primary 
CMB anisotropics, which play the role of white noise for 
the present purpose. 

In Fig. [31 we plot directly the mutipoles Cf~^ for the 
same two models. The effect of neutrino masses is clearly 
visible for all Z > 2 at Zm = 0.1, while for z^ > 1 it is 
necessary to reach Z > 20 in order to see a difference 
(since the maximum free-streaming scale is seen under 
a smaller angle at higher redshift). Remembering that 
the effect of neutrino masses on large I's can be split 
in two contributions, a matter power suppression and an 
excess of ISW, it is clear from the previous discussion that 
the latter effect contributes at all redshifts, but its most 
obvious manifestation is the fact that CT'-^ increases with 
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FIG. 3: Dimensionless cross-correlation spectrum in multipole 
space, l{l + V)Cj° for the same cosmological models as 
in Fig. [2] (i.e., with three neutrino species sharing the same 
mass vtiv = or niv — 0.41 eV). 



for large Vs. However, we will see in Sec. Ill Ci that only 
the region with I < 100 can be probed by observations: 
then, the neutrino- induced ISW effect is significant, but 
smaller that the opposite suppression effect. 

In Fig. 31 we plot the corresponding angular correla- 
tion functions 'uF'^{6). In this representation, the fine- 
structure of the high-^ multipole spectrum is by construc- 
tion averaged out, and it is not possible to see an ampli- 
fication at high Zm and small 6. The suppression caused 
by neutrino masses is visible for z„j = 0.1 at 6* < 15°, 
and for > 1 at 6* < 2°. 

In all these plots, we used only the linear perturba- 
tion theory. Doing so, the angular cross-correlation func- 
tions depend on the matter power spectrum inside the 
linear regime. To prove it, we compute again vo^'^{9) 
from the non-linear power spectrum obtained by apply- 
ing HALOFIT corrections [i^ to the linear one. The result, 
superimposed in Fig. HI is indistinguishable from that of 
linear theory. This shows that non-linear effects on the 
evolution of matter perturbations has much less impact 
than that of adding a neutrino mass. This is also true for 
the multipoles Cf^^ excepted for the smallest redshifts 
and highest Vs (for z,„ = 0.1, non-linear effects become 
important for / > 100). 



C. Detectability 

For a set of full-sky CMB and LSS experiments measur- 
ing the temperature multipoles aj^ (resp. galaxy-density 
multipoles af^) with a noise spectrum Nf (resp. iVp), 
the cross-correlation spectrum Cf^ can be reconstructed 
from the estimator 
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with a variance af'~^ given by 

. TG^2 (gp)^ + (gr+ivr)(gp^+^r) 

) = 2IT1 



FIG. 4: Angular cross-correlation function, multiplied by the 
average CMB temperature and displayed in units of micro- 
Kelvins, for the same cosmological models as in Figs. I2I3I 
We also plot the same functions including non-linear (NL) 
corrections to the matter power spectrum: they are indistin- 
guishable from the linear ones. 



Note that the estimator is not Gaussian, especially for 
small Vs: so, is only an estimate of the true (asym- 
metric) error bar on the reconstructed power spectrum. 
If the cross-correlation map can be reconstructed only 
inside a fraction /sky of the full sky, in first approxima- 
tion af^ should be multiplied by f^ky^'^- The variance is 



further reduced by v Al in case of data binning with bin 
width {Al). Note that in this case the covariance matrix 
is no longer diagonal, but nevertheless using a diagonal 
matrix under these approximations has been shown to 
work well, compared to the exact treatment, if we choose 
an adequately large binning (47j . In practice, for the 
multipole range in which we are interested, the CMB 
noise spectrum N'j'"'^ is much smaller than C'f" for ex- 
periments like WMAP and beyond, and can be safely 
neglected in the above expression. For a LSS survey con- 
sisting in a catalogue of discrete objects (galaxies, clus- 
ters, etc.), the noise spectrum is usually dominated by 
the shot noise contribution Nf^^ ~ 1/iV, where TV rep- 
resents the mean number of objects per steradian. The 
largest ongoing/future surveys (e.g. SDSS) should reach 
typically the order of 10® or even 10^. 

In Fig. O we show the typical errorbar that could 
be expected from a cross-correlation map with coverage 
/sky = 0.65 (corresponding to the usual galactic cut in 
CMB maps), using an ambitious LSS survey with surface 
density N = lO^st"^ in each redshift bin. We assumed 
b{z) ^ 1 for simplicity. These assumptions correspond 
essentially to the best measurement that could ever be 
done, since for such a high surface density the variance 
of the estimator of a single multipole product af^a'j^-^ is 
not affected by instrumental noise, and reduces to 
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(13) 



This expression can be interpreted as the product 
of the cosmic variance term Cf"^ times an enhance- 
ment factor depending on the correlation coefficient 



6 



{Cf^fKCl^C?^)- At large Vs, the late ISW contribu- 
tion to the total temperature anisotropy becomes vanish- 
ingly small, and the primary anisotropy plays the role of a 
large noise term, which cannot be removed. In this limit, 
the correlation coefficient is much smaller than one, and 
the variance aj^ gets correspondingly enhanced. Fig. [5] 
shows that the spectrum Cf^ can be reconstructed to 
some extent only in the range I < 100; beyond, one could 
only derive upper bounds. Note that the errorbar for 
each bin is roughly of the same order of magnitude as 
the effect of neutrino masses when f^, varies from to 
0.1. In Fig. [5l we also show the error degradation when 
/sky is reduced to 0.25 and iV to 7 x lO^st"^ in each red- 
shift bin. Finally, in Fig. [51 we plot the corresponding 
error bars for w'^^ (6) . Note that the synthetic error bars 
for w^^{d) are correlated with each other, unlike those 
for Cf^'^ . On small angular scales 6 < 1 (where the effect 
of neutrino masses is maximal) the Icr error on w'^'^{9) 
is of the order of 25%. 

We conclude from these estimates that the 
temperature-galaxy correlation power spectrum Cf^ 
is potentially sensitive to the neutrino mass in the 
observable range 10 < Z < 100, as well as the angular 
correlation function w'^'^{9) for 6* < 5° at z = 0.5 or 

< 3° at z = 1. Unfortunately, the enhancement of the 
ISW effect due to the impact of massive neutrinos on 
the linear growth factor is not directly visible: it would 
require precise data at high / and high redshift, for which 
the late ISW effect is masked by primordial anisotropics. 
The net effect of massive neutrinos on the observable 
part of Cf^ and on ViF^{9) is a suppression, caused by 
the usual free-streaming effect. However this effect is 
non-trivial in the sense that Cp*^ and Cf^ depend on 
f,y through different relations, due to the fact that the 
ISW term involves a time-derivative of the gravitational 
potential while the galaxy overdensity does not. Hence, 
the galaxy-temperature correlation spectrum can bring 
some information on neutrino masses which is not 
already contained in the sole galaxy auto-correlation 
spectrum. In the next section, we will quantify this 
statement by performing a parameter extraction from 
mock data accounting for future experiments. 



III. AN MCMC ANALYSIS OF MOCK DATA 



For a given data set consisting in various maps (i.e. 
multipoles a^) covering a fraction /sky of the full sky 
and assumed to obey Gaussian statistics, the likelihood 
function L is often approximated as 



/:cxH,i (detC;"')-^/2exp 
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where C°^^ is the data covariance matrix defined by 
[C;°'='^]xy = (aL^L)' and Cf^ the assumed theoreti- 
cal covariance matrix for a given fit, which contains the 



sum of each theoretical power spectrum C^^^^ and of 
the instrumental noise power spectra , estimated 

by modeling the experiment. Of course, the data co- 
variance matrix reconstructed from the observed maps is 
also composed of signal and noise contributions. Simu- 
lating a future experimental data set amounts in com- 
puting the noise spectra , given some instrumental 
specifications, and generating randomly some observed 
spectra C^^^"^^, given the theoretical spectra C/^^^'^ of 
the assumed fiducial model and the noise spectra . 
However, for the purpose of error forecast, it is sufficient 
to replace simply C^^ohs ^.j^g g^^^^ C/^^fid + N^^: 
this just amounts in averaging over many possible mock 
data sets for the same model, and does not change the 
reconstructed error on model parameters [48j . 

For instance, if one wants to estimate future errors for a 
CMB experiment, the maps to consider are temperature 
and ii^-polarization: X S {T, E} (here, for simplicity, we 
consider models with no gravitational waves and discard 
B-polarization) . The covariance matrices then read 



^obs 



th 



TBth 



(jTEt\\ 

^eeL^^ee 



(16) 
(17) 



Should one consider the combination of CMB data 
with a future galaxy redshift survey decomposed in N 
maps associated to N redshift bins, the matrices would 
become 2 + N dimensional, with an extra block 



I 



1, 



,N, (18) 



as well as non-diagonal coefficients [C;]i^2+i = C™' ac- 
counting for the late ISW effect. Note that all non- 
diagonal coefficients have no noise term, since the noise 
contributions in two different maps are expected to be 
statistically uncorrelated at least at first order. 

Finally, the option which is most interesting in our 
context, is to assume that the galaxy density auto- 
correlation maps arc not known (or just not considered, 
because they could be plagued by some systematic ef- 
fects), and that CMB data are only combined with the 
cross-correlation data, i.e. with N observed power spec- 
tra spectra ^ This is exactly what is being done 
in the current literature, in which authors try to get some 
new independent bounds on J^a from CMB plus CMB- 
LSS cross-correlation data, without employing LSS auto- 
correlation maps. In the approximation of Gaussian- 
distributed C; ' with central value C; ' and covari- 
ance given by 



(15) [Gov,],, ^ ((C™' - (Cr^')) (C™^ - (C™^))) 
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GiGjth , 



{21 + l)/sky 



(19) 
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FIG. 5: 68% error forecast on the power spectrum Cf^, which is displayed for the same two cosmological models as in Figures[2] 
(with = or 0.1). The smallest error boxes assumes a LSS survey with sl^y coverage /sky = 0.65 and surface density 
iV = lO^st"^ in each redshift bin. The binning in multipole space can be read from the width of each box. The largest error 
boxes correspond to /skyO.25 and N = 7 x 10** st^^. 
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FIG. 6; 68% error forecast on the angular correlation func- 
tion w'^^ , which is displayed for the same two cosmological 
models as in Figures [2j |3] (with — or 0.1). The error 
bars assume a LSS survey with sky coverage /sky = 0.65 and 
surface density N = lO^st"^ in each redshift bin. The spacing 
between each error bar reflects the binning width chosen in 
angular space. 



the likelihood of the cross-correlation data reads 



C(xUi (det Covi)-i/2exp 



I [Cov,]^iAf 



(20) 



with Aj = C™'°^"-C;^ The total likelihood is then 
the product of the CMB and cross-correlation likelihoods. 

In this section, we will focus on three ambitious future 
experiments: the Planck satellite, to be launched in 2008, 
which is expected to make the ultimate measurement 
of CMB temperature anisotropics, dominated by cosmic 
variance rather than noise up to very high I; the Dark 
Energy Survey (DES); and the Large Synoptic Survey 
Telescope (LSST), designed primarily for a tomographic 
study of cosmic shear, which would provide as a byprod- 
uct a very deep and wide galaxy redshift survey (close to 
ideal for the purpose of measuring the CMB-LSS cross- 
correlation since 

^G,G, ^ (jG^G, jgg^g^ Jqj. jnultipoles 
I < 100). For Planck, we computed the noise- noise spec- 
tra, Nf^ and iVf ^, hke in Ref. jH], with nine frequency 
channels. For the DES-like survey, we followed Ref. 



-iTGith 



and assumed a total number of galaxies of 250 million in 
a 5000 square degree area on the sky (or /sky = 0.13), 
with an approximate 1-cr error of 0.1 in photometric red- 
shifts, divided in four redhsift bins with mean rcdshifts 
Zi G |0.3, 0.6, 1, 1.3|, with the same selection functions 
as m Ref. m- For LSST, we used the same modeling 
as in (sH . with a net galaxy angular number density of 
80 per square arcminute and a coverage of /sky = 0.65. 
The galaxies are divided into six redshift bins with mean 
redshifts z, e {0.49,1.14,1.93,2.74,3.54,4.35}. For each 
bin the selection function, estimated bias bi and galaxy 
density are provided in [sij (Fig. 2, Eq. (16) and Ta- 
ble I) . The noise spectra iV; ' ' are then simply given by 
l/n^. 

We used the pubhc code COSMOMC [H] to do a Monte- 
Carlo Markov Chain (MCMC) analysis, fitting the theo- 
retical galaxy-temperature correlation to the mock data. 
For this purpose, we have written a module which com- 
putes the correlation multipoles following Eq. @ and the 
likelihood of the mock data given each model as described 
above. 

We then ran our modified version of CosmoMC for 
a model with eight parameters: the usual six of mini- 
mal ACDM (baryon density ^Ibh^, dark matter density 
^dmh^ , angular diameter of the sound horizon at last 
scattering 6, optical depth to reionization t, primordial 
spectral index ti^, primordial amplitude log[10^*'As]) plus 
the total neutrino mass rrii, and the equation-of-state pa- 
rameter w. Our fiducial model was close to the WMAP 
best- fitting model with = and w ~ — \. We consid- 
ered three possible combinations of data: Planck alone, 
Planck plus its cross-correlation with DES or LSST (but 
no information on galaxy auto-correlations), and finally 
Planck plus LSST, using all information and including 
the correlation. The probability of each parameter is 
displayed in Fig. [7] for each of these four cases called re- 
spectively CMB (Planck), CMB+GT (Planck+DES or 
Planck-hLSST) and CMB+GT+GG (Planck+LSST). 

Obviously the combination CMB+GT+GG does a 
much better job than CMB+GT for constraining all pa- 
rameters (and most spectacularly w and m^). This is 



0.0216 0.022 0.0224 

a if 




0.09 0.1 0.11 




1.039 1.0398 1.0406 




0.94 0.95 0.96 




2.96 2.98 3,„ 3.02 3.04 

log[10"AJ 





0.4 0.6 O.E 
£2 



0.5 1 1.5 

m (eV) 



FIG. 7: Marginalized probability of cosmological parameters 
obtained by fitting some mock data mimicking the properties 
of Planck, DES and LSST. The solid black curves accounts for 
CMB only, the red dashed for CMB+GT with Planck+DES, 
the blue dashed for CMB+GT with Planck+LSST, and the 
dotted blue for CMB+GT+GG with Planck+LSST (these 
combinations are precisely defined in the text). In each of 
these cases the cosmological model consists in ACDM (six 
parameter) plus an arbitrary total neutrino mass and 
equation-of-state parameter w. So, only eight of the above 
nine parameters are independent (as a consequence the prior 
for Ha is non-flat). The mock data is based on a fiducial 
model with rrii, — and w — —1. 
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FIG. 8: Two-dimensional marginalized likelihood contours 
involving obtained by fitting some mock data mimick- 
ing the properties of Planck, DES and LSST. The solid 
black curves accounts for CMB only, the red dashed for 
CMB+GT with Planck+DES, the blue dashed for CMB+GT 
with Planck+LSST, and the dotted blue for CMB+GT+GG 
with Planck+LSST (these combinations are precisely defined 
in the text). For each case, the two lines represent the 68% 
and 95% confidence levels. 



mainly due to the fact that the GT cross-correlation 
is partly screened by primary temperature anisotropics, 
while the GG signal does not have such an intrinsic noise 
contribution. We even try to repeat the CMB+GT+GG 
analysis with all C;*^^' correlations set to zero, and found 
no noticeable difference, showing that most sensitivity 
comes from GG rather than GT terms. However, the 
comparison between CMB alone and CMB+GT is still 
interesting per se. In fact, we are dealing here with 
an idealized situation, but in the future the GG auto- 
correlation signal could appear to be plagued by various 
systematic effects. In this case, independent information 
coming from the cross-correlation signal alone might be 
a useful piece of evidence in favor of the preferred model. 
Also, if the galaxy bias turns out to be very difficult to 
estimate with high enough accuracy, one may adopt the 
point of view of using the GG signal to measure bias, and 
the CMB+GT signal to estimate the best-fit parameters 
in some iterative scheme. 

In this prospective, it is interesting to note that the 
CMB+GT combination from Planck and LSST increases 
significantly the sensitivity of Planck alone mainly for 



^dmh'^ (by 30%), w (by 83%) and (by 38%). As a con- 
sequence, the sensitivity to the related parameter f^A in- 
creases by 76%. For our fiducial model with rrii, = 0, the 
95% confidence level upper bound on the total neutrino 
mass shrinks from 0.77 eV to 0.54 eV (for another fidu- 
cial model with 771^, > the sensitivity can only be larger 
than that, see e.g. ^^). At this level of sensitivity, the 
parameter is not correlated with Q\ or w, as can be 
checked by looking at two-dimensional marginalized like- 
lihood contours in Figure [H We conclude that the cross- 
correlation signal derived from Planck and LSST would 
have some useful sensitivity to both neutrino masses and 
dark energy parameters. Instead, the correlation between 
Planck and DES does not bring significant new informa- 
tion with respect to Planck alone. 

Ichikawa and Takahashi performed a similar fore- 
cast for Planck and LSST (with slightly different spec- 
ifications), using a Fisher matrix analysis rather than 
MCMC approach. They find a smaller sensitivity of the 
cross-correlation data to neutrino mass than we do, pos- 
sibly because of the various approximations entering into 
the Fisher matrix approach. 
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IV. CONCLUSIONS 

We have studied here the possibihty to use the cross- 
correlation between CMB and galaxy density maps as a 
tool for constraining the neutrino mass. On one hand 
massive neutrinos reduce the cross-correlation spectrum 
because their free-streaming slows down structure forma- 
tion; on the other hand, they enhance it because of the 
behavior of the linear growth in presence of massive neu- 
trinos. Using both analytic approximations and numeri- 
cal computations, we showed that in the observable range 
of scales and redshifts, the first effect dominates, but the 
second one is not negligible. Hence the cross-correlation 
between CMB and LSS maps could bring some inde- 
pendent information on neutrino masses. We performed 
an error forecast analysis by fitting some mock data in- 
spired from the Planck satellite. Dark Energy Survey 
(DES) and Large Synoptic Survey Telescope (LSST). For 
Planck and LSST, the inclusion of the cross-correlation 
data increases the sensitivity to rrii, by 38%, w by 83% 
and ridmh'^ by 30% with respect to the CMB data alone. 
With the fiducial model employed in this analysis (based 
on eight free parameters) the standard deviation for the 
neutrino mass is equal to 0.38 eV for Planck alone and 
0.27 eV for Planck plus cross-correlation data. This is far 
from being as spectacular as the sensitivity expected from 
the measurement of the auto-correlation power spectrum 
of future galaxy/cluster redshift surveys or cosmic shear 
experiments, for which the predicted standard deviation 
is closer to the level of 0.02 eV, leading to a 2a detection 
even in the case of the minimal mass scenario allowed 
by current data on neutrino oscillations (see [s^ for a 
review). However, the method proposed here is indepen- 
dent and affected by different systematics. So, it remains 
potentially interesting, but only if the neutrino mass is 
not much smaller than rrii, ~ 0.2 eV. 
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Appendix: the Limber approximation 

Let us consider some maps X{n) expanded in spherical 
harmonics 

oc / 

with 

at = J d^nY;,^{h)X{n). (22) 

The two-point correlation function of any two statisti- 
cally isotropic quantities X and Y can be expressed in 
terms of the power spectrum in multipole space 

Cr-iaLal:). (23) 

or in therms of the angular correlation function in a Lcg- 
endre polynomial basis {pi) 

^xr(^)^^2r^^^(^^^^)^xy, (24) 
I 

In the frame of observations, a direction dependent quan- 
tity X{n) is usually a quantity integrated over the line of 
sight, X{n) = J drX{x). The expression for a^, ((22|) . 
can then easily be transformed to Fourier space. Subse- 
quently expanding the plain wave in spherical harmonics 
and applying the completeness relation for spherical har- 
monics, one arrives at 

dr^X{kMkr)Y,l{k), (25) 

where X{k) is the Fourier transform of X{x), ji{r) is 
the spherical Bessel function, and k ~ \k\. This ex- 
pression can be simplified using Limbers approximation, 

/ dxf{x)ji{x) ~ y^2lTT / dxf{x)S{l ^ - x), leading to 
aL - (-)'/^ / 'i'-^X {k,k) YLikl (26) 
separating the k dependence of X in k and k — — — . 
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